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Abstract 

We formulate the basic theoretical methods for Bose-Einstein condensation of 
atoms close to a Feshbach resonance, in which the tunable scattering length of the 
atoms is described using a system of coupled atom and molecule fields. These in- 
clude the Thomas-Fermi description of the condensate profile, the c-field equations, 
and the Bogoliubov-de Gennes equations, and the Bogoliubov excitation spectrum 
for a homogenous condensed system. We apply this formalism to the special case 
of Bragg scattering from a uniform condensate, and find that for moderate and large 
scattering lengths, there is a dramatic difference in the shift of the peak of the Bragg 
spectra, compared to that based on a structureless atom model. The result is com- 
patible with the experimental results of [JJ. 

1 Introduction 

A recent experiment on Bragg scattering of a Bose-Einstein condensate (BEC) close to 
a Feshbach resonance Q], investigated the effect of varying the interatomic scattering 
length on the position of the Bragg spectral peak. This experiment used values of scatter- 
ing length, much larger than have previously been used in Bragg scattering experiments, 
and noted that there was no theoretically justifiable method of describing a BEC in this 
regime. 

In Paper I |2| we formulated a method for treating ultracold atoms close to a Fesh- 
bach resonance as a system of coupled atom and molecule fields. This method is ap- 
plicable in the regime covered by the experiment, and has been developed in such a a 
way as to be directly applicable to a c-field description of a BEC in the vicinity of a Fes- 
hbach resonance. The aim of this paper is to formulate the standard results, such as 
the Thomas-Fermi profiles and the Bogoliubov spectrum, appropriate to this model of 
systems of condensed atoms and molecules, and describe the differences which arise 
compared to the corresponding theories of structureless atoms with the same scattering 
length. 

This paper is a precursor to our next paper, in which full c-field computations of 
Bragg scattering will be performed, in a way which is directly comparable with the ex- 
periment of [T|. Here, we apply this model to the case of ideal Bragg scattering from 
a uniform condensate, and compare this with the corresponding theory for a BEC of 



structureless atoms, as done previously in (3). By doing this we obtain a benchmark 
comparison of the behaviour of a Feschbach molecule treatment with that of the struc- 
tureless atom models, avoiding the numerous complexities which necessarily arise in 
experiments. Thus, we can isolate the effects which arise only from the dynamics of the 
molecular basis of the Feshbach resonance, from those which are forced on us by the 
practicalities of experimental procedure. 

2 Formalism 

The formalism used here is outlined in Paper I [2i| and is based on the c-field methods 
described in detail in 0]. 

2.1 Phenomenological Hamiltonian for atom-molecule system 

As introduced in Paper I, the Hamiltonian of a coupled atom-molecule system is given 
by 

H= J dxjv> t (x)[-^-+V a (x))v>(x) + t (x)[--^-+V m (x) + c)0(x) + 

^?-$\x)$\x)${x)${x) + I [$Hx)f(.x)f{x) + fHx)fHx)^>(.x))Y (1) 

where U a a - 4nh 2 a bg l m is the background interaction strength and V a and V m are the 
external trapping potential for the atoms and molecules respectively. The terms with 
factors Uam and U mm correspond to atom-molecule and molecule-molecule scattering. 
However, in the systems considered here, the molecule field arises only during collisions, 
and is very small, making these terms negligible. 



2.1.1 Values of parameters 

The parameters g and £ are the coupling strength and detuning in the Feshbach reso- 
nance respectively. In our formalism they are given in terms of the experimentally mea- 
surable binding energy and s-wave scattering length by 

h 2 a 2 (7i-2Aa s )[l-2Aa bg t(^)/7i) 
2m Aa s {l+ t{j))-7i 

8nh 4 a 2 [a bg [n - 2Aa s ) -na s ){l-2Aa bg t(j i )lii) 
m 2 2Aa s [l+ t(f ))-n 

where t{x) = x-arctanl/x, the parameter A is the momentum space cutoff, and h 2 a 2 /m 
is the molecule binding energy for a certain value of the s-wave scattering length a s [2)- 
It is important to note that in this phenomenological model the parameters £ and g take 
on different values, depending on the closeness to the Feshbach resonance. 

In this paper we will consider a BEC of 85 Rb, for which there is experimental data of 
the binding energy close to the Feshbach resonance at 155 G \5\. For a given momentum 
cutoff A we can then calculate the values of the coupling g, and detuning e, for each 
value of the scattering length, using equations (3) and (2). We fit curves to the values 
calculated from the experimental data in order to extrapolate to other values of the scat- 
tering length. The data, obtained using |2) and (3) are well fitted by a linear relationship 
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Figure 1: The Feshbach resonance detuning e, and coupling parameter g, as functions of the in- 
verse scattering length. The black circles indicate the values calculated using the experimental data 
from 1 5j and equations (3) and (2). The black lines are fits to these calculated values. The shaded 
blue areas indicate the range of scattering lengths of interest in the experiment of |T| , as well as in 
the later sections of this paper. The momentum cutoff is here chosen to be A = 10 6 m _1 . 



between g and a s , and a quadratic relationship between e and a s . We show the data 
and fits in Fig.Q] 



2.1.2 C-field equations 



Since the atom field is usually much larger than the molecule field in the situations we 
shall consider, we can set U a m and U mm equal to zero and the c-field equations of mo- 
tion corresponding to the Hamiltonian Q~| become 



ih— — = — - — + V a {x) \y/(x)+U aa \y/(x)\ 2 y/(x) + gy/(x)*(p(x), (4) 
at \ 2m j 

_dd>(x) j h 2 V 2 \ g , 

iti^r- = \-- + £+^WP+yW. (5) 

at \ Am j 2 

In practice, when solving these equations numerically, it is necessary to explicitly use a 
projector in order to restrict the wavefunctions to the momentum subspace below the 
cutoff A, as described in Paper I and implemented in Paper III. 

The value of the momentum space cutoff A will, in the case of numerical simulation, 
arise from the simulation grid; any such computation is restricted to a finite number of 
momentum space modes. It is also necessary to introduce a projector in order to avoid 
effects from aliasing [4) . Furthermore, a momentum cutoff is necessary in order for the 
pseudopotential approximation to be valid [6j. However, as we show in Paper I, the actual 
choice of value for the cutoff required for our simulations has only a very small effect on 
the values of the phenomenological Hamiltonian parameters. 



2.2 Thomas-Fermi Solutions 



Stationary solutions to the equations of motion (4j and (5) can be obtained by letting the 
time evolution of the wavefunctions be i[/{x, t) = t//(x) exp{i[i a t/h) and cp{x, t) = (p{x) exp{ip 
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Figure 2: Thomas-Fermi profiles for a spherically symmetric condensate of 85 Rb for two of the 
different scattering lengths reported in Ref. 0: a s = 856ao (left panel) and a s = 2133ao (right 
panel). The trapping frequency is oj x = 27i x 17.5 Hz, and the total particle number — counting 
each molecule as two atoms — is 10 5 . The blue solid lines are the atom Thomas-Fermi profiles 
from (8), the green solid lines are the molecule profiles given by (9), and the dashed black lines are 
the Thomas-Fermi profiles obtained from the Gross-Pitaevskii equation. The parameter xo is the 
length scale associated with the trap, given by xq = \/h/2ma> x . 



It is clear that \i m - 2\i a in order for the coupling terms to be time-independent, so we 
get 



2pi a (p{x) 



2vt2 



+ V a [x) y/{x) + Uaa\f{x)\ fix) + gf{x) (p{x), 



2m 

h 2 v 2 \ g 7 

■- — + £ + y m (x)Uw + W. 

4m 2 



(6) 
(7) 



Taking the Thomas-Fermi limit and solving for y/ and (p gives the Thomas-Fermi solu- 
tions 



\y s {x)\ 

(p s {x) 



LL a -V a {x) 



Uaa + g 2 l2{2n a -E-V m {x))' 
g 



2(2ii a -e-V m (x)) 



y 2 {x). 



(8) 
(9) 



Fig.[2]shows the Thomas-Fermi wavefunctions for a spherically symmetric condensate of 
10 5 atoms of 85 Rb. The Thomas-Fermi solution for the atom wavefunction is similar to 
that for the structureless model for the same atom number and scattering length (also 
shown in the figure), with the important difference that the denominator has a spatial 
dependence. As long as the scattering length is moderate (left panel of figure |2), the de- 
tuning e will be large and the molecular trapping potential will be negligible in compar- 
ison, making the Thomas-Fermi solution indistinguishable from the standard Thomas- 
Fermi solution. However, for very large scattering lengths (right panel of figure[2), the 
detuning £ will be small enough for the molecular trapping potential to be significant, 
making the atom-molecule Thomas-Fermi profile different from that of the structureless 
atom model. 

Also, note that in (8), since the potential V m {x) is normally negligible compared to 
2[i a - e, the Thomas-Fermi solution for the atomic field is essentially of the same form 
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Figure 3: Density for a spherically symmetric condensate of 85 Rb for two of the different scattering 
lengths reported in Ref. \5\: a s = 856flo (left panel) and a s = 2133«o (right panel). The system 
is the same as that in figure [2] For this system, the portion of atoms in the form of molecules is 
approximately 0.6% for a s = 856flo and 4.5% for a s = 2133ao- (Note: At first glance, it appears that 
the solution y/GPE nas fewer atoms, than given by the atom-molecule model, but, as can be seen 
in the right hand figure, the discrepancy in the centre is compensated by at the wings, where the 
extra factor of inr 2 becomes more significant.) 



as that for the GPE at the same chemical potential. However, for the same total num- 
ber of atoms — counting each molecule as two atoms — the chemical potential for the 
GPE is slightly different from that found in this model. As well as this, the molecule field 
corresponds to the number of elementary molecules, and each physical molecule is a su- 
perposition of an elementary molecule and an atom pair, as discussed in Paper I. Closer 
to the Feshbach resonance, the proportion of atom pairs can become more than 50%, 
although in the systems we study here, these effects will be so small that they can be 
neglected. 

The atom density of the condensate in the Thomas-Fermi limit is given by 

n{x) = \y s {x)\ 2 + 2\<p s (x)\ 2 . (10) 

Fig-El shows the density profile for a spherically symmetric condensate of 10 5 atoms of 
85 Rb in the Thomas-Fermi limit. For comparison we have also plotted the Thomas-Fermi 
density profile obtained from a structureless atom model for the same scattering length 
and atom number. 



3 Quasiparticle excitations 

In the Bogoliubov approximation, the field operator can be separated into a mean-field 
ground state and an operator part describing the excitations. In our coupled system, we 
expand the atom and molecule wavefunctions in terms of quasiparticle bases as 



y/(x,t) ) 




' yroix) ) 




u p (x) ) 


b p e-™v' + 


' v* p {x) ) 




K (f>{X,t) j 




y (po(x) j 




K r p {x) j 









(11) 



where the first terms are the condensate atom and molecule ground state wavefunctions, 
and bp and bl are the quasiparticle destruction and creation operators respectively. Un- 
like the standard Bogoliubov expansion, where the only mixture is between the atom 
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creation and destruction terms, in this case it is also necessary to include the molecule 
creation and destruction terms. The quasiparticle operators in this expansion thus have 
both atom and molecule components. 



3.1 Two-component Bogoliubov-de Gennes equations 

Making these substitutions to the equations of motion (4) and (5) gives, after collecting 
terms with the same phase, the two-component Bogoliubov-de Gennes equations 

LU p + MV p = hojpUp, (12) 
LV p + M f Up = -hcjpVp, (13) 

where U p - (Up(x),r p {x)) T and V p - {v p [x), s p [x)) T . Here L is Hermitian, but M need 
not be Hermitian. In our case they are given by 



L = 



^- + V a (x)-ll a + 2UaaVo(x) 2 gVoW 

-^L + £+V m {x)-2 l i a 



UaaVoix) 2 + g</>oM 

M = 





(14) 



(15) 



Both L and M are in fact Hermitian. When y/Q and (po are chosen to be real, they are also 
symmetric. The multi- component equations fT2l[T3t are similar to the single-component 
ones. Indeed, if we let g — «■ in these equations, the upper components obey the stan- 
dard Bogoliubov-de Gennes equation for a single atom field. 



3.2 Orthogonality and normalisation conditions 

The normalisation and orthogonality conditions of the quasiparticle amplitudes can be 
derived by using the symmetry properties of L and M, for details see the appendix. We 
get the conditions 



>p,p 



(16) 



J dx^Up-V^Vp) = S 
f dx(v*U p >-Uj;v p ,) = 0. (17) 

We can also use the Bose commutation relations for the components of y/ = {y/,(p) T , 

[VaM.vjjt*')] = S atP 8{x-x'), (18) 
[y/ a {x),y/ p {x')] = 0, (19) 

from which we get the conditions 

L[ u *.P u i P -Kp v P,p} = <W 
p 

L(Ua,P V lp-K,pUp, P ) = 0. (21) 
P 

The four conditions (T6] [T7] [20] [2TJ are related; indeed (20) and (2D can be derived from 
QD and Q7). 
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4 Uniform condensate 



If considering a uniform condensate, the Thomas-Fermi solutions (8] |9) are exact and 
give 



fa = 
(po = 



U aa + g z /2{2Ha-£) 

2 _ 



(22) 
(23) 



where we have chosen fo to be real and positive. 

For the uniform condensate, the quasiparticle amplitudes can be expressed as plane 
waves, according to 



ric,l(x) 



Ukje 



ikx 



n/V 



71* O 

u k,r 



rk.ie 



Vv 

ikx 



s/V 



r * „-ikx 
r k,r 



where V is the volume of the system. 

In this case we can express L and M as 



L = 



1 r ^F" fJ-a + 2U aa n a gs/rTa 



gV n a 



(24) 

(25) 
(26) 
(27) 



(28) 



M = 



IU aa n a + g^fn^ 0' 
0, 



(29) 



Since there will be two distinct eigenfrequencies for each value of the momentum k — 
one corresponding to a state which is mainly atomic, and one corresponding to a state 
which is mainly molecular, as we shall see in the following section — the subscript p 
refers to the different momentum modes k, as well as the two eigenvalues, which we 
refer to as the "atomic" cja and "molecular" %. This terminology corresponds to the 
behaviour of the eigenvectors for sufficiently small |fc|. The situation becomes somewhat 
complicated at higher values of |/c|, as we shall show in the next section. 



4.1 Eigenvalues of the Bogoliubov-de Gennes equations 

The eigenvalues of the Bogoliubov-de Gennes equations are now given explicitly by a 
somewhat intricate procedure as follows: 
A : Define the quantities 



3? = L 2 U - M u + £22- 

§ = (£f 1 -Mf 1 - ( C| 2 ) 2 + 4Cf 2 ((,Cn + /;22) 2 -Mf 1 ). 



(30) 
(31) 
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B : define the quantities K l and K as the smallest and largest positive values of | fc| for 
which S = 0. 

C : For the range of values of | k\ given by < | k\ < K , 



-^-yji-VE, (32) 

V2 



~\piWE. (33) 
V2 



D : For the range of values of | fc| given by K l < | fc| < K h , 



o>k,A = -=}/$.- iy/l§\, (34) 
v2 



w fc , M = — ^V^+ iVi§]. 

v2 

E : For the range of values of |fc| given by K h < |fc , 



0>k,A = -^V^+^> (36) 
v2 



0)k,M = -^=y0lW§. (37) 

Fig.Hshows the eigenvalues calculated using this procedure for a uniform condensate of 
85 Rb with a total density of n = n a + 2n m = 10 20 m~ 3 at a scattering length of 900 a®. Here 
we use the experimental values of the binding energy taken from (5! to determine e and 
g, using Equations (2) and |3j . The momentum is measured in units of the momentum 
q of the Bragg pulse in the experiment of Ref. [1] with q = 1.6 x 10 7 m _1 . The following 
points should be noted: 

1. As can be seen clearly in the figure, the atomic eigenvalue t>k,A is similar to that 
obtained from the Gross-Pitaevskii equation, but is slightly shifted from this, due 
to the coupling between the atoms and the molecules. This shift is not dramatic, 
but is still measurable, for example by using Bragg scattering, as in the experiment 

of m. 

2. The atomic eigenvalue ci>k,A belongs to the atom- dominated state. The energy to 
create an atom- dominated quasiparticle is thus ftw^. 

3. Similarly, the energy to create an molecule-dominated quasiparticle is given by 
h(i)jc : M- Since (i>k,M is related to the binding energy of the molecules, and £ is nega- 
tive, (n)k,M is also negative. Thus, the energy to create a molecule-dominated quasi- 
particle is negative, and increases with |fe|, which is expected for a bound state. 

4. The higher eigenvalue (i>k,M nas the same general behaviour as (i)k,A> but since it is 
related to the binding energy of the molecules, it is shifted downwards. As k —- 0, 
and as long as the density is moderate, (i>k,M can be approximated by 

ha) 0iH x£ + 2Uaan-4iJ, a , (38) 

and as n — it will approach the value of the molecule field detuning e. 
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Momentum k [q] 

Figure 4: Eigenvalues W/t.A (blue solid line) and &)fc,M (green solid line) as functions of the mo- 
mentum k for 85 Rb at 900flo- The atom- dominated eigenvalue is similar to that obtained from the 
Gross-Pitaevskii equation (grey solid line), whereas the molecule-dominated eigenvalue is related 
the molecule field detuning e (grey dash-dotted line). The total density of the condensate is here 
10 20 m" 3 and q is 1.6 x 10 7 m" 1 . 



4.1.1 High momentum instability 

As can be seen in Fig.|D there is a crossover which occurs when a>k,A + ^k,M - 0> where 
the energy to create simultaneously a molecular- dominated quasiparticle and an atom- 
dominated quasiparticle is zero. As one proceeds through the crossover, the eigenvalue 
formula is determined successively by the procedures C, D and E, as given above. The 
midpoint of the crossover region is at the momentum 

h 2 k 2 2 

— - = 2fi a --(2U aa n + e). (39) 
2m 3 

Fig.[5]shows the behaviour of the high and low eigenvalues in the crossover region, where 
each eigenvalue has an imaginary part. The maximum amplitude of the imaginary parts 
of the eigenvalues at approximately 1.5 kHz. For this system, we would therefore expect 
a relatively fast instability to occur — on a time scale of approximately 0.07 ms — with 
atom-molecule pairs being created and destroyed. 

The instability arises because the energy to create a molecular quasiparticle of mo- 
mentum k becomes equal to the energy to destroy an atomic quasiparticle of momen- 
tum -k. This is a natural instability to expect, and provides a mechanism for the con- 
densate, which is metastable when the underlying interaction is attractive, to achieve its 
true ground state, a condensate of molecules. The relatively high value of the momentum 
required would make this a very much less important phenomenon in a trapped conden- 
sate, since the wavefunctions of high energy quasiparticles are largely located outside of 
the condensate. Since the transformation can only take place where there is a conden- 
sate, this would significantly decrease the size of the imaginary part. 
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Figure 5: High eigenvalue (green dashed line) and low eigenvalue (solid blue line) for the same 
system as in figure [4] in the region where the energy associated with the destruction of an atom 
coincides with that of the creation of a molecule. The eigenvalues become complex where the 
overlap exists, with their real parts shown in the left panel and their respective imaginary compo- 
nents are shown in the right panel. As in figuref?] the grey line is the eigenvalue obtained from a 
system of structureless atoms. 



4.2 Eigenvectors of the Bogoliubov-de Gennes equations 



Using the normalisation condition 1201 we can find a unique set of eigenvectors of the 
Bogoliubov-de Gennes equations corresponding to the eigenvalues | |3ZH37) : 



Uk,A = 



u k,M 



1 A k,M^ l ~ C k,M^ B k,M 



2 

k,A 



k,M' 



k,M 



i2 p2 
k,M k,M' 



k,A' 



1 -^- a IaKa 
!-4m ' 



(40) 



(41) 



and, for I - L,H, 



Vk.A 
?k,A 
Sk,A 



= Aic t AUk,A, 

= B kiA u k}A , 

= ChAVjcA, 



(42) 
(43) 
(44) 



where 



A k ,A 

B k ,A 
Ck.A 



r 2 



ha) k>A -L22 
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f^k.A ~ &22 
£l2 
tw)k,A + £22 



(45) 

(46) 
(47) 



The components of the eigenvectors for the same system as that in Fig.|4] are plotted in 
Fig-El as a function of the scattering length, for low momentum (Fig. |6(a)) and moder- 
ate momentum (Fig. |6(b)| . At low scattering lengths, the eigenvector corresponding to 
the low eigenvalue a) k>A is dominated by the atom components u kiA and v k>A , with the 
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Figure 6: Eigenvector components u^a (blue solid line), v^a (green dashed line), a (red dash- 
dotted line) and sj- A (cyan dotted line) for the low eigenvalue Wjt.A Q ert P ane ') and the high eigen- 
value ft)j;,M (right panel). The momentum is measured in units of the momentum q of the Bragg 
pulse in the experiment of Ref. Q] with q= 1.6 x 10 7 m _1 . 



molecule components r^ ^ and Sk t A becoming more significant the larger the scattering 
length. For larger scattering lengths, the amplitude of the molecule destruction operator 
is clearly significant compared to the atom one, and there is a large portion of molecules 
in the atom state. 

For the higher eigenvalue ci>k,M' me situation is reversed with the dominating eigen- 
vector component being the one associated with the molecule destruction operator. Here 
we also see increased mixing between the two states as the scattering length is increased. 



5 Application to Bragg scattering 

The quasiparticle description of the excitations in our atom-molecule system can be 
used to describe the excitation spectrum from Bragg scattering. Following the proce- 
dure of ID and (7), we derive the momentum transferred onto the condensate from an 
Bragg pulse, and compare the results we obtain here to those in (3). 
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The equations of motion for the coupled atom and molecule fields subject to a Bragg 
pulse are given by 



ih— — = — - — + v a (x) \y/(x) + Uaa\y{x)ry/(x) + hv pty/ + gy(x)*(p(x), (48) 

at \ 2m ] 

^d(b{x) I h 2 V 2 \ , g 9 

in— = he+ V m (x) \Mx) + 2nV „ t cos{q- x- (Dt)(p+ — w (x), (49) 

dt \ Am I v 2 

where 

Vopt = Vo(t)cos(q-x-cot), (50) 



where q and w are the wavevector and the frequency of the Bragg pulse respectively 13 18 1 , 
and Vq is the amplitude of the optical potential, given in terms of the Rabi frequency O 
and the excited state detuning A, 

no. 2 

V = —. (51) 

The optical potential for the molecule is chosen to be twice that of the atom on the as- 
sumption that the atoms in the molecule are very weakly bound, and for these purposes 
behave almost independently. 



5.1 Quasiparticle evolution 

Using the expansion fi"D , and the Bogoliubov-de Gennes equations, we can write the 
equations of motion 

i^{u p b p e- la '' t + v;b* p e lw A = V opt yr, (52) 
P 

where 

(1 0\ 



"opt - ^opt 



V 2, 



(53) 

Similarly to the method in (5), we can project out b p from the above equation by writing 
f dx{[ul,V opt yf + V^,V opt yf*)}= (54) 
= f dx£{il/J, (Upbpe-™ p f + V^y"^) -iVj, {u* p b;e M P t + Vpbpe-^P 1 )}, 

(55) 

= iZ{jdx [u^Up - V^Vp) bpe-™? 1 + fdx [u^v; - V^U*] b* p ^, (56) 

= ib p ,{t)e~ m p'', (57) 

where on the last line we have used the normalisation and orthogonality conditions (T6] 
\T7] . We therefore have 

bp{t) = -le^P 1 J dx{[up 0VX y+vft m y*}}, (58) 
Setting y/ = {y/ a ,(po) T , we get 

6 p (t) = -ijT e^^'dt' J dx\v opt [y/ {u p +v* p ) + 2(p {r* p + s* p )]Y (59) 
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Assuming that the Bragg pulse is square, with amplitude Vq and duration T, we get 



b p {T) 



-iV e lw P TI2 J dx [y/ (u* + v*) + 2cp {r* p + s*)) 



-iwT72 \q 



sin ((top + to) 772) 



COr, - to 



ton + to 



(60) 



(61) 



5.2 Momentum transfer from uniform condensate 

If considering the ideal case of Bragg scattering from a uniform condensate, the quasi- 
particle amplitudes are given by |40]-|44), and the ground state is given by 02 [9). The 
expression for b p now simplifies to 

b p (T) = -^e la ^ TI2 J^dx(^{u;+v;) + 2^{r; + s;))e lkx 



-ko772 



sin ((«p- to) 772) 



+ e 



iwr/2 



sin ((top + to) 772) 



to,, - to 



= -iV (v^(«* + i;p+2V^(r*+5*)) 



£) i((Up-ai)r/2 



sin((top-w)r/2) 



+ e 



itftj p +ft))r/2 



sin (to p + to)r/2) 



(62) 



(63) 



where AT a and N m are the number of atoms and molecules, respectively. 

The total momentum imparted to the condensate can now be evaluated using the 
Bogoliubov expansions fTD , and the normalisation and orthogonality relations fl6l fTTI 

no, 



P{T) = ^ j dx^ifr*Vtfr+^(l>*V(p>+ c.c. 

= Y. h p\\ b p^\ 2 +\\ b p^f] 

< sin ((to - to) 772) \ 2 I sinCCto^/ + to) r/2) x 2 



3%lf 
2 / 



Wo) -to 



'TVoCw, + Vq) + 2\/N m {r q + s q )\ 



(64) 



It is clear that this expression will have its maximum values at the points where cj-Wqj, 
and if the atom field is much larger than the molecular field — as is the case for the 
system studied here — the only significant maximum will be at the lower eigenvalue (i>l- 
Consequently the Bragg resonance peak will occur at the eigenvalue of the atom state, 
d)k,A in figure |U corresponding to the momentum of the Bragg pulse. 

Fig.[7]shows the fraction of particles being scattered out of the condensate during the 
Bragg pulse for three different scattering lengths. The portion of scattered particles is 
related to the momentum transfer as P/hqN, where N= N a + 2N m is the total number of 
particles. The resonance peak of the Bragg spectrum shifts to higher frequencies as the 
scattering length is increased. However, this shift starts to decrease when the scattering 
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Figure 7: Bragg spectrum showing the fraction of particles scattered out of the condensate as a 
function of the Bragg frequency for a uniform system 85 Rb with total density of 10 20 m~ 3 . Here 
we show the spectra for three different scattering lengths: a s = lOOao (green dashed-dotted line), 
a s = 700flo (red dashed line) and a s = lOOOflo (blue solid line). The duration of the Bragg pulse is 
T = 0.1ms, and the intensity of the pulse is Vq = 2n x 1 kHz. 



length is larger than approximately 700<2o- The free particle resonance is located at l/2nx 
hq 2 /2m~15A kHz. 

Fig-El shows how the shift of the Bragg resonance from the free particle resonance 
changes with the increase in the scattering length. For comparison we have also plot- 
ted the result based on the model of a condensate of structureless atoms with increased 
scattering length |3]. The two models agree well for low scattering lengths, but differ 
dramatically for larger scattering lengths. In the figure we have also marked the scatter- 
ing lengths for which there is available data of the binding energy, as in |5]. The values 
outside of these are calculated by the method described in Section|2] 

6 Conclusion 

The Bragg peak shift obtained from our model of a coupled atom-molecule condensate is 
significantly different from that obtained by a model of structureless atoms. Even though 
it describes the ideal case of Bragg scattering from a uniform condensate and is there- 
fore not directly comparable to that obtained in the experiment of Ref. [TJ, it is still qual- 
itatively similar to this. Most importantly, unlike other approaches to this problem, for 
example that by Kinnunen et al 1101 . the lineshift in figure [8] shows a clear downward 
behaviour for large scattering lengths, exactly as was reported experimentally in [TJ. 

In Paper III we will implement the formalism described here and in Paper I in full nu- 
merical calculations of Bragg scattering from an inhomogeneous BEC. The results from 
these simulations will be directly comparable to the results from the experiment of (TJ . 
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Figure 8: The shift of the peak of the Bragg spectra from that of the non-interacting gas, plotted 
against the scattering length a s . The solid blue line is based on the peak position calculated by 
1641 , and the black dashed line is the equivalent calculation based on the model of structureless 
atoms, as in [3] . The black triangles mark the values of the scattering length for which Ref. [ 5 gives 
experimental data of the binding energy. 



A Orthogonality and normalization conditions 

The orthogonality conditions of the quasiparticle amplitudes can be derived by writing 



(65) 

ha> p , f dx [u^Up - vJVp) = fdx (u p ,L + V^M 1 ) U p + Jdx (vj£ + + t/J,M) V p . 

(66) 

The Hermitian nature of L means that the two right hand sides are equal so that, for an 
appropriate normalization we can write 





S 



(67) 



Since both £ and M are symmetric we can then write 





V p T Z + U T p M*)Up,-{ 



(68) 

t/ p r £ + VpM)v p ,}, 
(69) 



so that, assuming a) p ^ -a) p >, we can say 




(70) 
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We can also find normalisation conditions using the Bose commutation relations for 
the components of iff - {\f/,<p) T , 

8 a:P S(x-x') = \y/ a (.x),\ff^x')] (71) 

= £ {u a , p (x)U* p ,(x') [b p {x),tf p ,(x')\ + V* p {x)V Pip ,{x') [b p (x),b p ,{x')\ 

p.p' 

+ U a ,p(x)V Pipl (x') \b p {x),b pl {x')] + V* >p {x)U^ p ,[x') [hl(x),b p ,(x')]tyz) 



from which, assuming that the quasiparticle operators b p obey the usual Bose commu- 
tation relations, we get the normalisation condition of the amplitudes, 

T,( U <*,P U lp-K,pVp, P )=8 a , P > (73) 
P 

Similarly, the other commutation relation, 

= [y/ a {x),y/ p {x')] (74) 

= £ \u a ,p{x)Ufs tP >{x') [b p {x),b p ,{x')] + V* ip {x)V* p ,{x') [bl(x),b p ,(x')j 

p.p' 

+ U at p{x)V* p ,{x') \b p {x),b p ,{x')\ + Vl p {x)U Pipl {x') [P p {x),b p ,{x')^ (75) 
gives the condition 

EKp^ -Vl p U p J = 0. (76) 

P 
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